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' Using a Langevin description of spinodal decomposition in fluids, we examine domain 
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growth in the diffusive, viscous and inertial regimes. In the framework of this model, nu- 
merical results corroborate earlier theoretical predictions based on scaling arguments and 
dimensional analysis. 



in . 

0^ ■ The dynamics of phase transitions in binary fluids quenched into the coexistence region has been the 

^ ' subject of considerable study in recent years |^,^. It is generally accepted that long after the quench, the 

^ ^ ^ phase separation dynamics can be characterized by a single time dependent length scale, R(t) t"'. As a 

Q ' result, much attention has focused on how domains grow in time — specifically what is the growth exponent 
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' Scaling and dimensional analyses due to Siggia jsj, Furukawa San Miguel et al [|| and more recently 

Bray H address this question. Experimental |q| and numerical [R-Ddl studies, however, have not necessarily 



supported these theories, sometimes providing conflicting results |15[ |. Often overlooked in spinodal decom- 
position in binary fluids is that several stages of growth can occur, in each of which a different transport 
mechanism dominates. This fact has been remphasized in [p|,^ Ul^Jl7[ | . Individual experiments and numerical 
simulations typically access only a particular regime. Lacking has been a clear demonstration of 1) the 
existence of these distinct regimes within a single model and, subsequently, 2) quantitative results in these 
regimes which validate theoretical predictions |0-pi. In this Letter we address these points. 
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To simulate phase separation in a binary fluid, we used the Langevin model of Farrell and Vails [Q. The 
order parameter ip is the difference in the concentration of the two fluid components. Its evolution and that 
of the fluid velocity are given by 



where p is the average mass density, F is an order parameter diffusion coefficient, and rj is the shear viscosity. 
Here a = 77(1 — 2/d) + ^, where C is the bulk viscosity, and d is the spatial dimension. The dimensionless 
constant. A, couples the order parameter to the fluid velocity and is also the strength of the convective 
flow. The chemical potential, fi — where F is the free energy of the system at equilibrium given by 
u] — ^ J d'^r[pu'^ + ^aip'^ — bip"^ + PK\Vip\'^]. The strength of the interfacial energy is PK. Below the 
critical temperature, a and b are positive constants. 

After the fluid is quenched, single phase droplets form and grow. In the coarsening process a competition 
between hydrodynamic and thermodynamic effects can lead to three dynamical regimes: the diffusive, viscous 
and inertial We discuss these briefly, using dimensional analysis based on Equations (1) and (2). 

In the diffusive regime, the fluid velocities are small, and the advective term in is negligible compared 
to the order parameter diffusion. Therefore, (1) becomes dtip ~ T'V^p. Since the chemical potential p k/R 
where k is the surface tension and R is the characteristic length scale in the system (i.e. domain size), we 
have R(t) {TK)iti. The coefficient Tk implies that the growth in this regime is driven by diffusion and 
surface tension. In two dimensions, for example, the surface tension k is given by k — ^{2(3K)^/'^ Q. 

In the viscous regime, hydrodynamics becomes relevant. In particular, in the velocity equation (j2|), the 
viscous term dominates the inertial terms. If one ignores the inertial and bulk viscosity terms, the shear 
stress term is balanced by the force due to the gradient in the chemical potential. Thus, yyV^u ^ XipVp so 
that R{t) ^ ^t. This is the linear growth law predicted by Siggia The coefficient ^ indicates growth 
driven by the surface tension and controlled by the viscous force in the fluid. The length-scale, Rd, and 
time-scale, td, at which the system crosses over from the diffusion regime to the viscous regime is given by 



(1) 



pdtu = ryV^u + a\/{\7 • u) - AV • (puu) - XipVfj, 



(2) 



setting {TK)hl ~ ^td. Thus, td ^ and Rd ~ {^)^ . 
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In the inertial regime, inertial effects dominate over the viscous forces so that ^ Xip\7fi. This leads to 

-R(^) ~ (^)^^^i SIS predicted by Furukawa [Q (see also ||^). The coefficient here indicates that the growth 

is driven by the surface tension and controlled by the inertial effects. The crossover between the viscous 

1 — 

regime and the inertial regime thus occurs at length scale and time-scale t/^ where ^ {^)'^th ^'^ 
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that th ^ pxif^i s-iid Rh ^ 'x^- Similarly, the length, Ri, and time, ti for crossover from diffusion directly 
to inertial are given by Ri ^ (^"x^)^ ^* ^ 1^^ respectively. This would correspond to the inviscid flow 
case. 

To facilitate growth of domains in each of these regimes and to access each of them within the framework 
of a single model, we vary Rd and Rh (relative to system size) by adjusting the parameters rj, A and p. For 
convenience, the actual (dimensionless) numerical equations we solve are the following: 

dt(l> = V2[03 _ ^ _ ^y2^J _ . J^^J ^ 

^tv^ = fi\/^v, + a^V.VfcWfc- A0V,;[(/>3 - (/.-/? V20] - A^[Vfc(w,Wfe) +VfcV,Ufc] + Wi. (4) 

k k 

The rescaled order parameter and transport quantities are given in terms of those used in (1) and (2) by 
(t> = (f)^, V = {■^)iu, 7] = (7 = A = A(p5^)^. Space and time are rescaled by r ^ r, 

1 -2 

t TKt. The dimensionless crossover lengths are given by Rd = (f ^^'^ ~ f^' ^^^'^^ have set 
p = 1. In two dimensions, the surface tension, k — |(2/3)2, so that by varying (3, we can control the surface 
tension. 

We studied deep, critical or symmetric quenches with (0) = throughout the course of the simulations, 
where (.) denotes an ensemble average or space average. The order parameter and velocity are initially 
taken as Gaussian fields with {(p) = (u,) = 0, and (0^) — (vf) — 0.005. The grid size Ax used was 1.7 
and the time step At was chosen as 0.05 in two and 0.02 in three dimensions, respectively. The mmrerical 
integration scheme is the same as in 0,0. The average domain size was defined as the first zero of the 
equal time correlation function G(r, t) = (0(x, t)(j){x. + r, <)), the Fourier transform of which is the structure 
factor, S{k,t). The fields, /i and Wi, were Gaussian, white noise with covariance given by the fiuctuation- 
dissipation relation [7|JTl|]. We found that adding noise does not alter the growth exponent in the scaling 
regime. However, it introduces curvature in the early growth so that longer times are required to reach 
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the scaling regime. The results we report here were obtained in the absence of noise and in all cases were 
averaged over 3 or 4 independent runs. 

In two dimensions, on a 1024^ system, we let ?) = 1, A = 1 and (3—1. Thus, Rd ^ Rh ^ ^ (in lattice 
units) are both small compared to the lattice size, L = 1024, so that for domain size R{t) ^ Rh, the system 
will favor droplet growth in the inertial regime. The data represented by (□) in Fig.l shows that R{t) has 
behavior consistent with a — 2/3. In order to have a viscous regime, one requires Rd <^ R{t) ^ Rh- This 
is satisfied by choosing, for example, r) — 20, A = 1 and (3 — 1 so that Rh ^ 120 and Rd ^ 3. In Fig.l 
the symbol (x) shows the growth under these conditions. It is consistent with a = 1/2 growth over a time 
interval spanning about 1.5 decades. The exponent of a = 1/2 in two dimensions was predicted in Since 
Rh ^ 120, the inertial force would not be expected to influence the growth until at late times when R{t) is 
comparable with Rh- To indicate how the a = 1/2 growth could yield to the a = 2/3 growth, we changed 
parameters to make Rh smaller so that the crossover from the viscous regime to the inertial regime can 
happen earlier. The symbols (O) and (+) in Figure 1 show data for Rh ^ 30 {fj — 11, X = 1, (3 = 1) and 
Rh ^ 7 iv — 5, X — 1, f3 — 1) respectively. As Rh (and Rd) decreases, the data shows that the viscous 
growth and a later, faster inertial growth occurs progressively earlier. Finite size effects and the need for very 
long times to see adequate viscous and inertial growth make quantitative analysis of growth in the crossover 
regimes difficult. The time evolution of the Reynolds number. Re, the ratio of inertial to viscous effects, is 
consistent with the behaviour of R{t) as Rh decreases. Corresponding to the parameters for domain growth, 
the insert to Figure 1 shows how Re changes from its behavior in the viscous regime, where Re < 1 and is 
essentially constant, to that in the inertial regime where it increases as t^/^ @. We find that for 77 > 11, the 
system lies well within the viscous regime until the influence of inertial flow at very late times. The scaled 
correlation functions G(^) and C^G(^) § are shown in Figure 2, for fj — 20, X — 1, f3 = 1, where ^ = r{k) 

r }^i^(^f^ t)dk 

and (fc) = — - — ^- — . The data collapses well for several times, indicating that the a = 1/2 growth is in the 
scaling regime. The scaling behavior for a = 2/3 was shown in [ pl| . 

Using this model in two dimensions, we previously [ pT| examined the behavior of a as a function of the 
coupling constant, A(0 < A < 1), by fixing fj, [3 and p. For this one parameter system, Rh ~ ^ and 
Rd ~ TiT^- For 1/2 < A < 1, Rh ^ Rd ^ 1, so that only the inertial growth survived. For A ^ 0, the domain 



size R(t) < Rd and the dominant mechanism was diffusion. It was the first attempt to show within a singie 
modei different regimes, Iiowever, the one-parameter system was iimited in how weil it could capture all 
three regimes. The competing mechanisms of viscosity, inertial force and surface tension appear to demand 
a system with two parameters. Two-dimensional lattice Boltzmann and lattice gas simulations seem to be 
carried out primarily with relatively small Rh, thus the a = 2/3 estimates are consistent with growth in the 
inertial regime P, p^ , p^ . The results from molecular dynamics are controversial. It has been pointed out 



that the a — 1/2 growth obtained in |14| may be attributed to droplet coalescence. Velasco and Toxvaerd 
Ip^ observed a ~ 1/2 crossing over to a = 2/3 in their two-dimensional molecular dynamics simulations. 

Three-dimensional simulations were carried out on Equations (3) and (4) using a system with 256"^ lattice 
sites and show behavior analagous to that observed in two dimensions. As above, we set (3 — 1. If ?} = 1 and 
A = 1, Rh ^ Rd ~ 1 and are small compared to the domain size R(t). One thus expects inertial growth with 
a — 2/3 at late times, and this is seen by the data represented by (□) in Fig. 3. If f) = A (with Rd ^ 1), 
the system should favor growth in the viscous regime for sufficiently large fj. The symbols (x) and (O) in 
Fig. 3 show growth for tj — 25, A = 25, and ij — 20, A = 20, respectively. As fj increases, the growth becomes 
consistent with a = 1. The crossover between the viscous regime and the inertial regime can be simulated 
through decreasing Rh, while keeping Rd small (~ 1). The symbols (+) in Fig. 3 used ?} = 12, A = 12 
and show that a regime with a growth exponent of 1 gradually yields to a slower growth regime, a 2/3 type 
growth. Figure 3 (insert) shows that the behaviour of the Reynolds number. Re, is consistent with growth 
for appropriate parameters in the inertial and viscous regimes. Finite size effects are more pronounced in 
three dimensions, so that the inertial regime is difficult to access as Rh increases. In Fig. 4 is plotted the 
scaling of the correlation functions G{^) and ^^G'(^) (insert) for the inertial regime in three dimensions. The 
quality of the collapse of the data in Fig. 4 for several times indicates that the a = 2/3 growth is in the 
scaling regime. 

Earlier work by Farrell and Vails on the same model was carried out on an Sl'^ lattice with t) 1, A ^ 1, 
and a — 2 so that Rh ^ 1. Their estimate of a 1 was based on an extrapolation of a time-dependent, 
effective exponent in terms of inverse droplet size. Puri and Dunweg used a Cell Dynamical System 
model and obtained a ^ 1 on a model (with 80'^ lattices) without the convective term in the velocity equation 
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and with ry = 1, A ~ 2 ,and a — 1. Using their parameter values with the convective term on a 128^ lattice, 
we find an early a 1 growth that crosses over to a slower a 2/3 growth at later times. Shinozaki and 
Oono [|j and Koga and Kawasaki obtained a ~ 1 at late times with their models (Model H), ignoring 
the inertial terms. It was noted in |^ that for larger values of viscosity there is a crossover from a ~ 1/3 
to a ~ 1 growth. Such a crossover can occur because a larger viscosity increases Rd which then favors 
the diffusive growth for domain sizes R[t) < Rd- Lattice Boltzmann simulations also provide linear growth 
estimates |9|Jl0|. The model we have used allows for slight compressibility. However, it has been shown that 
in the viscous regime a does not change with the incompressible condition To our knowledge, the 

a = 2/3 in the inertial regime has not been observed in experiment or three-dimensional simulations. 

In summary, we have used a single model system to probe the hydrodynamic regimes that a phase sepa- 
rating fluid can undergo. In particular, we have shown how domain growth can be favored to take place in 
these regimes by an appropriate choice of the crossover lengths Rd and R^ within a finite size simulation. 
Moreover, we have obtained values for the growth exponent a in these regimes in two and three dimensions 
that are in agreement with the predictions of scaling and dimensional arguments. Our work helps to explain 
the estimates of growth exponents, a, obtained in a number of previous studies. 
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support from the National Science and Engineering Research Council of Canada. Y.Wu is grateful to the 
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using the computational resources of the Advanced Computing Laboratory at the Los Alamos National 
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Figure Captions 

FIG. 1. The domain growth R{t) vs time t in two dimensions for different crossover lengths Rh showing 
the change from growth in the inertial regime (t^/^) to the viscous regime (t^^^)- The errors in the 
data due to different initial conditions are of the order of the size of the symbols. The insert shows the 
Reynolds number, Re — vR{t)/fi, as a function of time, where v is the characteristic velocity calculated 
as dR/dt Q. The data are consistent with Re ^ t^/'^ in the inertial regime and Re ~ constant in the 
viscous regime. The straight line has slope 1/3. The symbols represent fj = 20, A = l,Rh ^ 120(x), 
?) = 11, A = 1, Rh 30(0), 77 = 5, A = 1, i?,, - 7(+) and r) = 2, A = 1, Rh l(n). The surface tension 
controlling parameter /3 = 1 for all cases. 

FIG. 2. The scaled and normalised correlation function G(^) vs ^, where f = r{k) and the first moment 

r kSik t)dk 

(fc) = •' n , for ?7 20, A = 1, and /3 = 1 on 1024 x 1024 lattices, at t = 2,500 (O), 3000 (x), 

4000 (o), 5,000 (+) and 6,000 (□). The insert shows f^G(^) versus ^ for the same times and since 
(k) 25, scaling is good for ^ 80 lattice units. 

FIG. 3. The domain growth R{t) vs time t in three dimensions for values of crossover lengths Rh showing 
the change in a from inertial {t'^^^) to viscous (t) regime. The errors are of the size of the symbols. 
The insert shows that the behaviour of the Reynolds number. Re, is consistent with t^^^ in the inertial 

regime and t in the viscous regime. The straight lines have slopes 1/3 (- - -) and 1 ( ). The symbols 

represent = 25, A = 25, (o), 7^ = 20, A = 20, (O), = 8, A = 8, (+), ?? = 1, A = 1, (□). The surface 
tension controlling parameter /3 — I for all cases. 

FIG. 4. The scaled and normalised correlation function G(^) vs f for ry = A = 1 and /3 = 1 on 256^ lattices, 
at t = 600 (o), 800 (□), 1,000 (x), and 1,200 (o). The insert shows C^G(0 vs ^ for the same times. 
Since (k) ^ 6, the times scale well to ~ 40 lattice units. 



